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The variation of the period distribution function of late-type binaries is studied. It is shown that the Taurus-Auriga pre-main sequence 
population and the main sequence G dwarf sample do not stem from the same parent period distribution with better than 95 per cent 
confidence probability. The Lupus, Upper Scorpius A and Taurus-Auriga populations are shown to be compatible with being drawn 
from the same initial period function (IPF), which is inconsistent with the main sequence data. Two possible IPF forms are used 
to find parent distributions to various permutations of the available data which include Upper Scorpius B (UScB), Chameleon and 
Orion Nebula Cluster pre-main sequence samples. All the pre-main sequence samples studied here are consistent with the hypothesis 
that there exists a universal IPF which is modified through binary-star disruption if it forms in an embedded star cluster leading to a 
general decline of the observed period function with increasing period. The pre-main sequence data admit a log-normal IPF similar to 
that aiTived at by Duquennoy & Mayor (1991) for main sequence stars, provided the binary fraction among pre-main sequence stars is 
significantly higher. But, for consistency with proto-stellar data, the possibly universal IPF ought to be flat in log-P or log-semi-major 
axis and must be similar to the Kl IPF form derived through inverse dynamical population synthesis, which has been shown to lead 
to the main sequence period function if most stars form in typical embedded clusters. 
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1. Introduction 



The initial distribution of orbital periods of binary-star systems, 
i.e. the initial period function (IFF), poses an important con- 
straint on star-formation theory as well as being a necessary in- 
put for modelling stellar populations. Constraining its form, and 
its possible variation with star-forming conditions, is thus of fun- 
damental importance. Observations of late spectral-type systems 
(stellar masses < 2 Mq) have shown that many young popula- 
tions have a higher binary proportion than old Galactic-field pop- 
ulations (Mathieu 1994; Zinnecker& Mathieu 2001 for reviews; 
and more recently Duchene et al. 2007a, 2007b, Connelley et 
al. 2008), at least within a certain range of binary star periods. 
Duchene (1999) has performed a comprehensive comparison of 
various young stellar groups in terms of their binary properties 
and concluded that very young star forming regions are very 
likely to have a binary excess compared to the field or star clus- 
ters. This conclusion is enhanced by the study of proto-stars by 
Connelley et al. (2008). 

However, variations between young populations also appear 
to be evident, such as comparing the Taurus-Auriga (Kohler & 
Leinert 1998) and Orion Trapezium Cluster populations (Prosser 
et al. 1994, Petr et al. 1998, Petr 1998, Simon et al. 1999) that 
have a similar age (» 1 Myr). Connelly et al. (2008) even study 
the possible variation of the shape of the proto-stellar period 
function as proto-stars transcend the various stages of their spec- 
tral energy distribution evolution. The source of such variations 
must be understood in order to unearth any possible systematic 
variations with the boundary conditions of their birth environ- 
ment. 



In this contribution we postulate the existence of a universal 
IPF and test it against observations of binary period distributions 
for various star forming regions. 

1.1. Binary Formation 

The formation of binary systems remains an essentially unsolved 
problem theoretically (e.g. Maury et al. 2010). Fisher (2004) 
shows analytically that isolated turbulent cloud cores can pro- 
duce a (unknown) fraction of binary systems with the very wide 
range of orbital periods as observed. But direct cloud collapse 
calculations are very limited in predicting binary-star proper- 
ties owing to the severe computational difficulties of treating the 
magneto-hydrodynamics together with correct radiation transfer 
and evolving atomic and molecular opacities during collapse. 
Available results applying necessary computational simplifica- 
tions suggest the preference of a typical period, arouncQ IP ^ 4, 
but such a narrow period distribution cannot be transformed 
to the observed wide range of periods, Q < IP < 10 (Kroupa & 
Burkert 2001). 

The currently most advanced hydrodynamical simulations 
have been reported by Moeckel & Bate (2010). They allow a 
turbulent SPH cloud to collapse forming a cluster of 1253 stars 
and brown dwarfs amounting to 191 Mq. The cluster has a radius 
of about 0.05 pc and contains a substantial binary and higher- 
order multiple stellar population with a large spread of semi- 
major axis, fl, but peaking at a few AU. After dynamical evolu- 
tion with or without expulsion of the residual gas the distribution 
of orbits ends up still being quite strongly peaked at a few AU 
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with a significant deficit of orbits with a > lOAU when com- 
pared to the main-sequence population (their fig. 11). This state- 
of-the art computation therewith confirms the above stated is- 
sue that it remains a significant challenge for star-formation the- 
ory to account for the Gaussian-type distribution of a spanning 
10"' - 10'* AU as observed by Duquennoy & Mayor (1991) and 
Raghavan et al. (2010) for G-dwarfs, by Mayor et al. (1992) for 
K-dwarfs and by Fischer & Marcy (1992) for M-dwarfs in the 
Galactic-field. One essential aspect which is still missing from 
such computations is the stellar feedback which starts heating 
the cloud as soon as the first proto stars appear. These heating 
sources are likely to counter the gravitational collapse such that 
in reality the extreme densities are not achieved allowing a much 
larger fraction of wide binaries to survive. 

More general theoretical considerations suggest that star- 
formation in dense clusters ought to have a tendency towards 
a lower binary proportion in warmer molecular clouds (i.e. in 
cluster-forming cores) because of the reduction of available 
phase-space for binary-star formation with increasing temper- 
ature (Durisen & Sterzik 1994, hereinafter DS). On the other 
hand, an enhanced hinnry proportion for orbital periods IP < 5.6 
may be expected in dense clusters due to the induction of binary 
formation through tidal shear (Horton, Bate & Bonnell 2001, 
hereinafter HBB), thus possibly compensating the DS effect. The 
initial period distribution function (IPF) may thus appear simi- 
lar in dense and sparse clusters, apart from deviations at long 
periods due to encounters and the cluster tidal field. A certain 
fraction of stars form in small-A^ systems, and the dynamical 
decay of these is likely to affect the final distribution of IP ^5 
binaries (Sterzik & Durisen 1998), giving rise to non-uniform 
jet activity (Reipurth 2000), but again, quantification is next- 
to-impossible given the neglect of the hydrodynamical compo- 
nent. But, the binary formation channel must be vastly domi- 
nating over the formation of non-hierarchical higher-order mul- 
tiples because otherwise even low-density pre-main sequence 
stellar populations would have a large fraction of single stars 
because initially non-hierarchical multiple systems decay on a 
core-crossing time scale 5x10'* yr, Goodwin & Kroupa2005). 
The low preponderance of high-order multiple proto-stars in ob- 
servational surveys is noteworthy in this context (Duchene et al. 
2007a). 

1.2. On the existence of a universal IPF 

Assuming binary systems emerge, by whatever process, with a 
universal IPF consistent with the pre-main sequence data, and 
that most stars form in embedded clusters of modest density, 
Kroupa (1995a, hereinafter Kl; 1995b, hereinafter K2) shows 
that the IPF evolves on a few cluster crossing time-scales to the 
observed Galactic-field, or present-day period function (PDPF). 
Applying this approach to dense embedded clusters, Kroupa, 
Petr & McCaughrean (1999, hereinafter KPM) and Kroupa, 
Aarseth & Hurley (2001, hereinafter KAH) find consistency with 
the observed period function in the Orion Nebula Cluster despite 
assuming the universal K I IPF. This suggests that observed vari- 
ations of the pre-main-sequence-star period function may be at- 
tributed to the dynamical history of the respective population, 
and that there may formally exist a universal parent IPF that 
emerges from the binary star formation process and from which 
the observed pre-main-sequence and main-sequence cases can 
be generated. 

Note that this universal IPF is strictly speaking not observ- 
able due to the rapid evolution after binary star birth, but the 
formal concept of a parent IPF is useful for synthesising binary 



star populations. This concept is equivalent to that of a universal 
IMF. Given that the parent distribution of stellar masses which 
emerges after the proto-stellar phase has been found to be invari- 
ant, it is natural to state the Invariant Initial Binary Population 
Hypothesis (Kroupa 2011). Indeed, given that it is known that 
the stellar IMF is invariant over a large region of physical star- 
forming regime it would be expected that a universal IPF ought 
to exist also. This is because the mass of a star is the last quan- 
tity to be established in the star-formation process, with the bi- 
nary property of a stellar system being established, essentially, 
at a prior formation stage. Discarding the notion of a universal 
IPF, the above star-cluster work at least shows that the dynami- 
cal history of a population needs to be taken into account before 
discussing variations of period distributions. 

There are two general processes driven by fundamentally dif- 
ferent physics that evolve an IPF with time. On the one hand, 
system-internal processes such as tidal circularisation, star-disk 
interactions and decay of small-A^ proto-systems affect binaries 
with periods IP ^5. K2 referred to this collectively as eigenevo- 
lution of the orbital parameters. Eigenevolution leads to period- 
eccentricity and period-mass-ratio correlations. On the other 
hand, binaries with IP>5 are affected by encounters in em- 
bedded clusters {stimulated evolution of the period function). 
Stimulated evolution does not significantly evolve the eccentric- 
ity distribution of a population, but leads to a depletion of the 
initial period and mass-ratio distribution function for compan- 
ions with weak binding energy (e.g. Kl; K2; Parker et al. 2009). 

The aim of this contribution is to quantify the significance 
of the variation of the observed period distribution of various 
populations of binary systems in order to allow an improved as- 
sessment of a possible underlying parent distribution of orbital 
periods or equivalently of semi-major axes. In doing so, the fol- 
lowing questions will be tested: Is there an invariant IPF? And, 
which form does this IPF have? Two IPF forms will be stud- 
ied, namely the Kl IPF (Eq.[4|below) which had been suggested 
previously to unify the pre-main-sequence and main-sequence 
binary populations, and a log-normal IPF. 

The IPF conceived here is a formal initial period distribution 
function established after the proto-stellar phase is finished. As 
such, it may not be observable, because eigenevolution aff'ects 
the shorter-period binaries on time scales of about 10'' yr, while 
stellar-dynamical and tidal-field influences affect it at larger- 
period on the dynamical time-scale of the star-forming system. 
The IPF is therefore as much a theoretical construct as the IMF 
is, which is also not observable (see the Cluster IMF Theorem, 
Kroupa 2008). The notion of a formally invariant IPF is very 
powerful. If this notion is verified to be consistent with obser- 
vational data it would be important as a boundary condition for 
star-formation theory (e.g. Fisher 2004), and also for the genera- 
tion of theoretical initial stellar populations such as for A^-body 
experiments of embedded clusters and 'classical' binary popu- 
lation synthesis (e.g. de Kool 1996; Johnston 1996), and for dy- 
namical population synthesis (K2; Marks & Kroupa 201 1). 

1.3. Overview 

Section |2] introduces the data sets collected from the literature 
and used in this contribution for the statistical analysis. Then, in 
Section|3]the significance of the difference between the observed 
canonical pre-main sequence and canonical main sequence pe- 
riod distribution is tested. In the next step of the present analy- 
sis (Section |4]i all data sets are tested for consistency with the 
Kl IPF and the log-normal period function of Duquennoy & 
Mayor (1991, hereinafter DM). All data sets are combined in 
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Section |5] to address the question if there exists a single parent 
period function of the Kl or the log-normal type, and a revised 
IPF is constrained using only 'matching' pre-main sequence data 
in Section|6l The conclusions follow in Section]?] 

2. Data 

Brief profiles of the observational data sets used here are pre- 
sented in the sub- sections below. In general, the binary star sur- 
veys provide for each detected binary system the observed pro- 
jected separation, A (in AU), which we transform to an estimate 
of the true semi-major axis, a, via a - A/0.95 (Leinert et al. 
1993). The values of a are then converted to an orbital period 
using Kepler's third law, P^^^ - a^/Mtot, where Pyis = P/365.25 
is in yr and a in AU, assuming the system mass Mtot = 1.3 Mq. 
This transformation is valid only statistically and not on an 
individual-object basis. Note that we choose a system mass of 
1 .3 Mq to reflect the primary being of solar mass with a typical 
secondary star which we take to be of about the average stellar 
mass. Leinert et al.(1993) used a system mass of 1 Mq, but this 
difference amounts to only about 10 per cent in log-mass and is 
thus of no consequence for the present analysis. In the case of 
the DM main-sequence sample and the pre-main-sequence bina- 
ries with IP <A (data by Mathieu 1994 and Richichi et al. 1994), 
the published IP values are adopted here. 

For each observational sample of A^obs.sys systems (each of 
which may be a binary or a single star) the thus obtained list 
of IP values is binned into nb bins over the interval IP\ to IP2 
covered by the observational survey (obtained from the survey 
interval in AU, a\ to a2, using Kepler's third law as above). The 
so constructed histogram has a bin-width 6lPi - (IP2 - lPi)/nb 
(in most cases equal for all bins except where noted otherwise), 
and per bin ; of width SIP, there are Nobs,bmilPi) binaries. The 
observational estimate of the true period distribution function, 
expressed as the binary proportion per unit IP interval, is thus 
obtained from 

fohsilPi) = A?obs,bin(/f ,)/(A^obs,sys SIP,), (1) 

with an associated binomial uncertainty (Petr et al. 1998), 

(foUlPi)[l-foUlPi)] ' 



eobs(/^,) = 



A^obs.sys SlPi 



(2) 



Binomial uncertainties are used here rather than Poisson uncer- 
tainties, ^eobsilPi) = (/obs(^f'/)/[A^obs,sys(^^^'/])'^^ since the typi- 
cal situation is at hand of observing NobsMnilP d successes in A^sys 
trials, each trial having a probability filPi) of success, whereby 
fobSPi) is an estimate of /(/P,). Note that eab^Pi) < ^e^bsilPi)- 

2.1. Main-sequence data (canonical; can. ms) 

The data by Duquennoy and Mayor (1992) is defined as the 
canonical main sequence sample because it comprises until re- 
cently the largest, and most detailed, multiplicity survey among 
low-mass main-sequence stars. It mainly consists of G-dwarf 
stars, but K- and M-dwarf main-sequence binaries have es- 
sentially the same period distribution (Fischer & Marcy 1992, 
Mayor et al. 1992, and fig. 1 in Kl). The characteristics of the 
DM study that qualify them for the 'canonical' main-sequence 
sample are the large and complete sample of solar-type stars 
within 22pc of the Sun, and the fact that all periods are covered. 
A key finding of the DM study is that the frequency distribution 
of the binary stars' orbital periods is unimodal and shows a well- 
defined log-normal distribution. In this analysis their published. 



incompleteness corrected, histogram of IP values, which is based 
on A^obs.sys - 164, is adopted. The range in period covered in the 
study is -1.0 < /P < 10 with SIP; = 1.0, and nb ^ II. 

A more recent survey of the multiplicity properties of solar- 
type stars has been presented by Raghavan et al. (2010). This 
survey comprises a sample which is about 2.5 times larger than 
the DM study and reaches to larger distances. In contrast to the 
DM survey, it is not based on a consistent 13-year spectroscopic 
survey but uses a compendium of observations from various 
sources. Raghavan et al. derive a period-distribution which they 
fit very well by a log-normal form with peak at logioP = 5.03 
and standard deviation of o-{logioP) = 2.28. The DM log-normal 
distribution peaks at logioP - 4.8 and has a standard deviation 
of cr(logioP) - 2.3. The Raghavan et al. study leads therefore to 
very similar results, and for the time being we continue to use 
the DM survey as the canonical one. 

2.2. Pre-main-sequence data (canonical, can.pms): 
Taurus-Auriga 

The pre-main sequence binary star data, i.e. the 'canonical' PMS 
sample, is composed of results obtained from multiple star sur- 
veys in the Taurus-Auriga star forming region. 

At short orbital periods, i.e. for IP < 2, binaries detected by 
the spectroscopic survey of Mathieu (1992, 1994) are consid- 
ered. Visual binaries in Tau-Aur, with 4.5 < IP < 7.5, have been 
extensively searched and studied by Ghez et al. (1993), Leinert 
et al. (1993), and Kohler & Leinert (1998), and the data from 
the latter two publications are included in the 'canonical' PMS 
sample. Slightly shorter period binaries (with 3 < IP < 4) have 
been addressed by Richichi et al. (1994), using lunar occulta- 
tion measurements. Hence, the 'canonical' PMS sample covers 
a significantly large range of orbital periods, from IP = 0.5 
to IP - 7.5, although not continuously (see lower panel in 
Figure. [TJ. All stars included in the 'canonical' PMS sample are 
low-mass stars, similar to the stars included in the 'canonical' 
MS sample. The PMS sample includes classical T Tauri stars and 
weak-line T Tauri stars that show indistinguishable binary star 
period distributions (Kohler & Leinert, 1998), and are thus con- 
sidered as one population in our 'canonical' PMS data. Table [T] 
lists the details for each dataset in terms of, e.g. how many tar- 
gets, Nohs.sys, have been observed, which are the lower and upper 
survey limits (expressed in lower and upper semi-major axis), or 
what distance to the region has been assumed. 

In the context of a possibly dynamical evolution of the bi- 
nary stars' period it is important to note that the average sur- 
face density of young stars in Taurus-Auriga is relatively low 
(Gomez et al. 1993). Variations or changes in the overall pe- 
riod distribution due to encounters over the age of the current 
Taurus-Auriga PMS population, which is ~lMyr, are expected 
to be small. Only during the earliest evolutionary phases of the 
systems, i.e. just after the formation of the stars, is some dynam- 
ical evolution conceivable (Kroupa & Bouvier 2003). 



2.3. Other Pre-main-sequence data 

Beside Taurus-Aurigae, the orbital period distribution of binary 
stars surveyed in the star forming regions of Lupus, Chamaeleon, 
Upper Scorpius, and the Orion Trapezium Cluster (Table [1] are 
considered here. By necessity the binary fractions per separa- 
tion range as published by the respective researchers are adopted 
such that here no re-analysis of possible contamination is per- 
formed as it would go beyond the scope of this contribution. 
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Table 1. Details of binary survey data used in this paper 



Star-forming region 


^ obs.sys 


nb 


Lower limit 


Upper limit 


Adopted distance 


References 








«i [AU] 


fl2 [AU] 


[pc] 




Tau-Aur (can.PMS) 


91 


2 


lp=0.5 


lp=1.5 


140 


(a) 


Tau-Aur (can.PMS) 


30 


1 


lp=3.5 


lp=3.5 


140 


(b) 


Tau-Aur (can.PMS) 


178 


4 


19.5 


1832 


140 


(c),(d) 


Lupu.s 


127 


4 


18.3 


1348 


160 


(e) 


Chamaeleon assoc. 


77 


4 


11.6 


460 


160 


(f) 


UScA 


68 


3 


18.0 


460 


145 


(g),(h) 


UScB 


51 


3 


18.0 


460 


145 


(g),(h) 


Orion TC 


62 


1 


214 


625 


450 


(i),a) 




109 


1 


72.9 


214 


450 


(i),a) 



Notes to table: Coloumn (2): no. of targets surveyed. Coloumn (3): no. of bins over the survey interval [IPi^IPt]- Column (4): lower limit of the 
survey interval in [AU], or alternatively in IP. Column (5): upper limit of the survey interval in [AU], or alternatively in IP. Column (6): distance to 
the region. Column (7): references; (a)= Mathieu 1994, (b)= Richichi et al. 1994, (c)=Leinert et al. 1993, (d)=Kohler & Leinert 1998, (e)=Kohler 
(priv. comm.), (f)=Kohler 2001, (g)=Brandner & Kohler 1998, (h)=Kohler et al. 2000, (i)=Petr et al. 1998, (j)=Petr 1998. 



Like for the canonical Tau-Aur PMS data, the surveys com- 
prise low-mass stars with a mixture of weak-line TTauri and 
classical TTauri stars. Although it cannot be ruled out com- 
pletely that some older ZAMS stars might be included in the 
samples, the large majority of the surveyed targets have been 
confirmed to be of PMS nature (e.g. Covino et al. 1997). 
Moreover, the orbital period distributions of the samples do not 
change significantly if the samples would be restricted to con- 
firmed PMS stars (e.g. Kohler 2001). None of the surveyed PMS 
populations in the analysed regions are older than a few Myr The 
Orion Trapezium Cluster, with an age of ~ 1 Myr (Hillenbrandt et 
al. 1997), is likely the youngest region, while Upper Scorpius is 
the oldest with ~5 Myr (Preibisch et al. 2002). The star forming 
regions Chamaeleon and Lupus have ages in between the Orion 
TC and Upper Scorpius [Chamaeleon: 2 Myr (Luhman 2004), 
Lupus: -2-5 Myr (Wichmann et al. 1997, Makarov 2007)]. 

The surveyed regions do show significantly different charac- 
teristics with respect to their average stellar densities. Nakajima 
et al. (1998) analysed the average surface density of compan- 
ions, taken as an indicator for the strength of clustering, and 
find that Lupus is the least clustered region, while the Orion TC 
is strongly clustered. The stellar density in the Orion TC was 
found to be ~ 5 X 10"* stars/pc^^ (McCaughrean & Stauffer 1994). 
Chamaeleon is considered as a loose stellar aggregate, similar 
to Lupus, but with a somewhat higher average stellar density 
(Nakajima et al. 1998). The region Upper Scorpius is part of the 
Sco-Cen OB association and does not show a very high stellar 
density like the Orion TC, although it is very likely that the orig- 
inal configuration of Upper Sco was much denser than today, but 
the stars have been dispersed with time. The stellar populations 
in Upper Scorpius are typically divided into Upper Scorpius A 
(UScA) and Upper Scorpius B (UScB), as these two regions 
show spatially distinct distributions. Both regions are found at 
very similar distances (de Zeeuw et al. 1999) and we assume the 
same overall distance value of 145 pc as adopted in the binary 
survey paper of Kohler et al. (2000). A notably difference be- 
tween UScA and UScB is that UScA contains several high-mass 
(B-type) stars, while no high-mass stars are present in UScB. 
Interestingly, the observed binary period distributions in UScA 
and UScB are different, with UScB showing a strong preference 
for wider binaries, while in UScA mostly binaries with small 
separations are present (Brandner & Kohler 1998). 

It should also be noted that Upper Scorpius and Orion are 
regions of low-mass and high-mass star formation, while Lupus 
and Chamaeleon are sites of low-mass star formation only. 



Table 2. Histogram of main-sequence, fms(lPi)^ and pre-main 
sequence, fpmsilPd, binary fractions and adopted uncertainties. 



IP, 






U.SP.) 




0.5 


0.024 


0.012 


0.035 


0.038 


1.5 


0.049 


0.017 


0.072 


0.038 


3.5 


0.079 


0.022 


0.200 


0.080 


5.0 


0.095 


0.024 


0.164 


0.021 


6.0 


0.082 


0.022 


0.156 


0.021 


7.0 


0.076 


0.022 


0.124 


0.020 



A survey of interest in the present context is that by Ratzka 
et al. (2005) and Simon et al. (1995) of p Oph which has a den- 
sity in between that of Taurus- Auriga and the ONC. Kohler et al. 
(2000) note that UScA and p Oph have very similar binary distri- 
bution functions. We will return to p Oph and other clusters such 
as Pleiades and Hyades in a more detailed investigation using 
Nbody computations. 

3. Main sequence versus pre-main sequence 

In this section the statistical significance of the difference be- 
tween the observed main-sequence and the pre-main sequence 
binary star period distributions is estimated. Focus is on the 
'canonical' main sequence and pre-main sequence distributions 
as defined in Section |2] The binary proportion per unit logioP 
interval is given by Eq. [I]/(/P,) = where Nbin(lPi) is the 

number of binaries in the ;th log-period interval 6lPi {- 1 here), 
and A^sys is the total number of systems in the survey, whereby 
each single star and each binary count as a system. In order to 
obtain matching data sets, the main-sequence histogram plotted 
in fig. 1 of KPM is Unearly interpolated to obtain three ZP, bins 
with unit width for /F, = 5.0,6.0,7.0 (i.e. covering the interval 
4.5 < IP < 7.5). The resulting data are listed in Table |2] 

The well-tried statistic is applied to test the null hypothe- 
sis that both observed distributions stem from the same underly- 
ing parent distribution, by calculating 

^{fUlPd-fpu..{lPi)f 

AT = y ^ -0 — - (3) 

eUlPd + el^dPi) 

to obtain — 15.7 with v - 6 degrees of freedom. Obtaining 
such a large value, or larger, has a significance probability 0.01 < 



Kroupa & Petr-Gotzens: The initial binary-star period function 



5 



V < 0.02, so that one can be less than 2 per cent confident that 
the null hypothesis holds true. Excluding the IP < 4 bins where 
the main-sequence and pre-main sequence data agree within the 
uncertainties, and concentrating instead only on the more recent 
Taurus-Auriga data (Kohler & Leinert 1998), = 13.2, v = 3 
with 0.001 < !P < 0.01. The null hypothesis can thus be rejected 
with approximately 99 per cent confidence. It may be concluded 
that the two observed distributions stem from different parent 
distributions. 

To obtain an additional assessment of the confidence 
in this result, the Wilcoxon-Signed-Rank (WSR) test (e.g. 
Bhattacharyya & Johnson 1977) provides a welcome alternative. 
This is a non-parametric test, making no assumptions about the 
form of the underlying populations, such as there being a well- 
defined mean and variance, as opposed to using the statis- 
tic which assumes such structure in the underlying distributions. 
The WSR test assesses the likelihood of observing a certain frac- 
tion of the data asymmetrically about a reference data set. To 
construct the WSR statistic, the differences fpmsilPi) - fmsUPd 
are ordered according to their absolute values. These are ranked, 
and the ranks associated with the positive differences are added 
to form the test statistic T'^. The statistic is symmetrical, that is, 
the same result is obtained by considering the negative differ- 
ences. For the data in Table |2] - 2\ with n - 6 data points. 
Obtaining such a large or larger has a significance proba- 
biUty P = 0.016 (table 10 in Bhattacharyya & Johnson 1977). 
The null hypothesis is thus only supported with a confidence of 
1.6 per cent, confirming the above conclusion. Using only the 
data with IP > 5 gives « = 3 differences, which is too small 
for this non-parametric test to allow significant conclusions, in 
contrast to the test used above, because the latter rehes on 
additional information about the populations. 

The Kolmogorov-Smirnov test cannot be applied here nor 
later in this paper, because the data sets cover different IP ranges, 
and /p was estimated by DM after applying incompleteness 
corrections, that is there exists no list of complete IP values from 
which a cumulative distribution can be generated. 



4. Standard period functions 

Having established with a high level of confidence that the 
canonical pre-main sequence and main sequence period distri- 
butions are significantly different, one can proceed to begin in- 
quiring as to how the parent distributions of the two data sets 
may be described. In the literature two analytic forms of possi- 
ble parent distributions for the two data sets have been used. 



4.1. Initial period functions (IPFs) and present-day period 
functions (PDPFs) 

For pre-main sequence binary systems Kl suggests an IFF, 



MIP) = /pms T] 



{IP - IP^in) 
6 + {lP-lP„,,nf' 



(4) 



using the canonical pre-main sequence and main sequence data 
as constraints. This is not in contradiction to the results of the 
previous section since one takes into account through exten- 
sive A^ody modelling that the main-sequence distribution re- 
sults from the pre-main sequence distribution if most stars form 
in modest (embedded) clusters. Kl found rj = 3.5 and 6 = 100 
give good fits to the pre-main sequence data, and to the main 



sequence data after passing through a typical star clustej^ as- 
suming the overall primordial binary proportion is /p^s = 1 , and 
setting the minimum orbital period to 1 d { IP min = 0). As seen in 
Fig.[T]note in particular that fyJJP) is essentially flat for IP > 4.5 
(fl >20 AU). Eq. ID with the above parameters is referred to as 
the Kl IFF, and IPFs constructed with different values for the 
parameters as of the K type. 

A more elaborate model based on the above Kl IFF, but in- 
cluding eigenevolution giving the correct correlation between 
eccentricity, period and mass-ratio for short-period binaries, is 
constructed by K2. This model has IPj^m = l,?? = 2.5,6 - 45 
and is required to re-produce and predict period- and mass- 
ratio-distribution functions and the distribution of orbits in the 
eccentricity-period diagramme for realistic stellar populations. 
Such details are however not required in the present treatment 
and wouldn't lead to different results. A useful feature of Eq. 
is that it can easily be converted to a period-generating function 
(Eq. llbinKl), 



IP{X) = /P™i„ + 



(5) 



where Xe[0, 1] is 

jl^-f^l(lP)dlP = /pn,s. 



a uniform random variate, and with 
jiP ./iviv" - J pms- The maximum allowed log-period, 

IPmax, follows from Eq.|5]with X = I. Eq.|5]allows efficient con- 
struction of a pre-main sequence population, such as for A^-body 
calculations of the evolution of embedded clusters (e.g. KPM; 
KAH). 

For the PDPF the Gaussian distribution in IP describes the 
observed period distribution of local G-dwarfs very well (DM), 



foMilP) - /ms K 



1 



1 , IP-avlP ^1 
2'- trip > 



(6) 



where crlP and avlP are, respectively, the variance in IP and the 
average-/P, and J^^ ' fouilP) dlP - /ms is enforced by adjust- 

X+CO 
^ foMilP) dip - Kfms- DM measured an overall 

binary proportion of /n,s = 0.58, crlP = 2.3 and avlP = 4.8. 
This is referred to as the DM PDPF, and to period functions con- 
structed according to Eq. |6]but with different values for the pa- 
rameters as of the DM type. A simple period-generating function 
cannot be written down; instead the Box-Muller method for gen- 
erating IPs is resorted to (e.g. Press et al. 1994). 



4.2. The tests 

The aim is to test with which confidence the IFF (the Kl- 
distribution, Eq. |4] with rj - 3.5, 6 - 100 and /p^s = 1) 
or the PDPF (the DM-distribution, Eq. |6] with crlP = 2.3, 
avlP = 4.8 and /ms - 0.58) can represent each data set of 
Section |2] For each observational sample the theoretical dis- 
tributions are converted to matching histograms by generating 
Mh.sys = 10* periods and binning these into the same ZP-bins 
as in the respective observational sample, obtaining (Eq. [1]) 
MlPi) = foN,h,binilPd/iNth.sysSlPi), where MlPd is either of 
the K (/o = /pms) or the DM (/, = /ms) type. Given the size of an 



^ The typical star cluster is the birth site of most stars in the Galaxy. 
According to Kl it contains typically 200 binaries in a characteristic 
radius of about 0.8 pc, while Adams & Myers (2001 ) find that most stars 
would originate from compact 10-100 member groups. This estimate 
can be argued to correspond to the previous one if residual gas loss 
and subsequent expansion with loss of stars from the modest embedded 
cluster are taken into account. 
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observational sample, A^obs.sys, the expected number of binaries 
in each IP bin is 



N,,p,bUlPi) = MlPi) A^obs.sys SlPi (7) 

with associated expected binomial uncertainty 

eN.exp.binC/i', ) = [MlPi) ( 1 - MlPd) 6lPi A^obs,sys] ' ■ (8) ^ 



To quantify the goodness-of-fit the statistic for the K or DM 
distribution becomes. 



nh I 



A^exp.bin(^f /) - Nob,Mn(lPi) 
^N,exp. 



(9) 



with nb degrees of freedom since no free parameters are fitted. 
Terms with NsxpMnilPd - give infinite since a finite da- 
tum is inconsistent with the model, which is treated as having no 
intrinsic uncertainties. The significance probability, !P, of obtain- 
ing Af^ > Xk('^^Xi)m^' i^ evaluated using the incomplete gamma 
function (Press et al. 1994). 

4.3. Results 

The above procedure is applied to all the data and the results 
are shown in Fig. [1] The conclusion with very high confidence 
(at the 99.5 per cent level or better) is that the canonical main- 
sequence data are not consistent with the Kl IPF (upper panel), 
and that the canonical pre-main sequence data are not consistent 
with the DM PDPF (lower panel). This does not change even if 
the IP - 9.5 main-sequence value is ignored in the upper panel 
(giving X- = 43 with P = 0.00 for the Kl IPF). This result is 
consistent with that of Section |2] where it was shown that the 
canonical main-sequence and pre-main-sequence data are sig- 
nificantly different. 

The results of applying this procedure to each of the pre- 
main sequence data sets are presented in Fig. |2] It is evident 
that the Taurus-Auriga, Lupus and UScA populations are con- 
sistent with the hypothesis that they be drawn from the Kl IPF. 
The Cham pre-main sequence population is marginally consis- 
tent with this hypothesis that the parent distribution is the Kl 
IPF at the 5 per cent confidence level, whereas the full UScB 
sample is only consistent with the Kl IPF at the 3 per cent con- 
fidence level. UScB becomes consistent at the 19 per cent confi- 
dence level with the Kl IPF if three of the widest binaries are re- 
moved from the sample, assuming these are significantly closer 
(< 90 pc) than the bulk of UScB stars at ~ 145 pc. Their appar- 
ent brightness is clearly brighter than expected for their spectral 
type when placed at 145 pc, and proper motion measurements 
confirm nearby distances for at least two of the wide binary sys- 
tems (Perryman et al. 1997; Salim & Gould 2003). 

On the other hand, an alternative hypothesis is that the data 
sets be drawn from the DM PDPF. This hypothesis can be dis- 
carded with very high confidence for the Taurus-Auriga, Lupus 
and UScB pre-main sequence populations. It cannot be dis- 
carded for the Cham, UScA and ONC populations. The Cham 
and UScA populations can, in fact, be drawn from either period 
function. The enhanced significance probability for the hypothe- 
sis that Cham be drawn from the DM PDPF nevertheless allows 
the hypothesis that Cham started off in a clustered mode with a 
Kl IPF, the Cham population possibly being an evolved version 
of the canonical pre-main sequence population. This also holds 
true for the ONC sample as shown by KPM and KAH. 
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Fig. 1. Testing the likelihood that the observed canonical binary- 
star period distributions fit the Kl IPF (dashed curve) or the DM 
PDPF (solid line). The model histograms are shown as dotted 
and solid lines with expected binomial uncertainties as error-bars 
(Eqs|2]and|8]l. Solid circles are the G-dwarf main sequence data 
in the upper panel and canonical pre-main sequence data in the 
lower panel (Section|2]i. The;i'- value and significance probabil- 
ity, f, of observing such a large or larger x^ is written in each 
panel (upper numbers for testing the data against the Kl IPF, 
lower numbers for testing against the DM PDPF.) 



5. Parent period distribution functions 

In Section|43]the Kl IPF and the DM PDPF were tested against 
all the available data sets individually. 

In this section data sets are combined and more general so- 
lutions for possible parent distributions of the K and DM types 
(Eqs|4]and|6]respectively) are sought by scanning the parameter 
spaces 77,(5 and crlP,avlP using suitable increments by evalu- 
ating Eq. |9] and the associated significance probability 'P. This 
procedure is demonstrated in Figs. |3]and|4]by first of all find- 
ing all admittable solutions to the canonical data sets. The solu- 
tions so found agree with the known solutions (the Kl IPF and 
DM PDPF), and one may be confident that the algorithm func- 
tions. 

On combining «dat data sets under a common hypothesis, the 
combineA x\g^y^ and number of degrees of freedom, n^comb, are 
computed. 
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Fig. 2. Similar to Fig. [T]but for all pre-main sequence samples. In all panels the dotted continuous curves are the model distributions (Kl IPF upper curves, DM PDPF lower 
curves), whereas the model histograms are shown as solid lines with expected binomial uncertainties as dotted error-bars. Thick histograms are the data (Section|2}. The;^^^ value 
and probability, P, of observing such a large or larger is written in each panel (upper numbers for testing the data against the Kl IFF, lower numbers for testing against the 
DM IPF). For UScB, the lower set of numbers refers to the dotted histogram, which results after removing the three widest binaries from the sample. 
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Fig. 3. Survey of rj, 6 parameter space for solutions of the K type 
with /pms = 1 (Eq. m using the canonical pre-main sequence 
data. Open circles and small dots delineate models consistent 
with the data at the 1 and 5 per cent confidence level or better, 
respectively, whereas the thick solid dot shows a model consis- 
tent at the 50 per cent level or better The Kl IPF is shown as the 
thick cross. Note that the Kl IPF is not the same as the best so- 
lution found here. This is not surprising because the Kl IPF was 
derived by constraining the IPF by both the pre-main sequence 
and main-sequence period distribution functions. The best so- 
lution found here thus lacks one major constraint which lead 
to the Kl IPF, but both are sufficiently close to conclude that 
the Kl IPF is a good solution to the pre-main sequence data 
alone. The straight line shows the asymptotic solution for large 
•5: fKiilP) = ir]/6) {IP - /P^n), with i]/6 = 0.0339. 



"dal 

/ 

where and nbj are the chi-square and number of degrees of 
freedom of data set /. Possible data sets are / - Lupus, Cham, 
UScA, UScB, ONC, canonical pre-main sequence (=can.pms), 
canonical main sequence (=can.ms)]. The significance proba- 
bility, 'pKorDM, of obtaining x~ ^ A'comb' evaluated as above 
(Eq. |9]), and 'Pk,'Pdm are the resulting combined confidence 
probabilities on using the K and DM type PFs, respectively, for 
testing the hypotheses set up in the following. A hypothesis is 
deemed consistent with the data if !PKorDM > 5 x 10"~, i.e. if the 
significance probability is 5 per cent or better 
The following hypotheses are tested; 

1. A PF of the K type is the parent distribution of all combined 
data sets («dat = 7). 

Assuming /p^s = 1 and scanning the parameter space, as 
above, yields no solutions at all in agreement with the data 
at the 0.001 per cent level or better (!Pk < 10"^). Thus the 
hypothesis can be discarded at the 99.999 per cent level. 
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Fig. 4. Survey of crlP, avlP parameter space for solutions of the 
DM type with f^s - 0.58 (Eq.|6]l using the canonical main se- 
quence data. Open circles and small dots delineate models con- 
sistent with the data at the 1 and 5 per cent confidence level or 
better, respectively, whereas thick solid dots show models con- 
sistent at the 50 per cent level or better The DM PDPF is shown 
as the thick cross. 



Assuming fp^s - 0.8 and f^s - 0.6 also leads to rejection of 
this hypothesis with 'P > 99.999 per cent confidence. 
The general conclusion, stated with a confidence better than 
99.999 per cent, is thus that there is no single parent distri- 
bution of the K type for all data sets combined. 

2. A PF of the DM type is the parent distribution of all com- 
bined data sets («(jat = 7). 

Assuming /^s = 0.58,0.8, 1.0 and scanning the parameter 
space, as above, yields no solutions at all in agreement with 
the data at the 0.001 per cent level or better (!Pdm < 10"^). 
Thus there is no single Gaussian parent distribution of all 
data sets combined. The hypothesis is rejected with a confi- 
dence better than 99.999 per cent. 

3. A PF of the K type is the parent distribution of the combined 
canonical pre-main sequence and canonical main sequence 
data sets (wdat = 2). 

The relevant confidence is determined from Eq. [10] with 
[/ =can.pms, can.ms]. The result is that there is no solution 
with confidence better than !Pk - 10"^ for/pms - 1,0.8,0.6, 
so that this hypothesis can be discarded with a confidence 
better than 99.999 per cent. 

4. A PF of the DM type is the parent distribution of the com- 
bined canonical pre-main sequence and canonical main se- 
quence data sets («dat — 2). 

The relevant confidence is determined from Eq. [TO] with 
[; =can.pms, can.ms]. There is no solution with confidence 
better than !Pdm = 10"^ for /^s = 0.6. For /„,,, = 0.8 the 
parameter region crlP - 2.5 + 0.4, av/P - 5.1 ± 0.2 has 
V < Q.l per cent, while for /ms = 1.0 the parameter region 
o-lP = 2.6 + 0.4, avlP = 5 .2 + 0.2 has !P < 0.05 per cent. This 
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Fig. 5. Testing hypothesis |5] with /p^s = 1-0. Survey of ri,S pa- 
rameter space for solutions of the K type with /p^s - 1 (Eq. |4]i 
to the combined canonical pre-main sequence data. Lupus and 
UScA pre-main sequence data sets. Open circles and small dots 
delineate models consistent with the data at the 1 and 5 per cent 
confidence level, respectively. The Kl IPF is shown as the thick 
cross, and the straight line is as in Fig. |3] The thick solid dots 
show models consistent at the 50 per cent level or better. 



hypothesis can thus be discarded with 99.9 per cent confi- 
dence. 

5. A PF of the K type is the simultaneous parent distribution of 
the canonical pre-main sequence data, the Lupus and UScA 
pre-main sequence data sets (wdat = 3). 

The relevant confidence is determined from Eq. [10] with 
[; =can.pms. Lupus, UScA]. These are chosen because Fig.|2] 
indicates these data sets to be similar in that they appear 
more or less unevolved according to the results of Kl and 
K2. Scanning rj, 6 parameter space leads to solutions with 
!Pk > 10--if/p„s = 1.0,0.8 (Fig. [5] and |6l respectively), 
whereas !Pk < 10"^ if /pms = 0.6. Thus, a simultaneous 
parent distribution is only possible if the pre-main sequence 
binary-star fraction is larger than the main sequence value 
of 60 per cent. 

6. A PF of the DM type is the simultaneous parent distribu- 
tion of the canonical pre-main sequence data, the Lupus and 
UScA pre-main sequence data sets («dat = 3). 

Scanning crlP, avlP parameter space leads to solutions with 
fuM > 10^2 if /ms = 1.0 (Fig. 13), whereas Pdm < 5 x lO^^ 
if /ms = 0.8 and Pdm < 10"^ if fms = 0.6. Thus, a simul- 
taneous Gaussian parent distribution is only possible, with a 
significance probability better than 5 per cent, if the pre-main 
sequence binary-star fraction is significantly larger than the 
main sequence value of 60 per cent (note that in Eq.|6]is 
here the binary-star fraction required for the combined pre- 
main sequence samples). 
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Fig. 6. As Fig.|5] Testing hypothesis[5]with /p^s - 0.8. 

7. A PF of the K type is consistent with all pre-main se- 
quence data (canonical. Lupus, UScA, Cham, UScB and 
ONC) (Hdat = 6). 

No parent distribution of the K-type can be found with 
0.1 per cent confidence or better This suggests that these 
pre-main sequence data sets may either have different IPFs, 
or that they stem from the same IPF but were modified, as 
would be expected for example for the ONC data (KPM, 
KAH). UScB, however, forms a definite outlier, since it has 
a significant surplus of binary stars at long periods (Fig. |2), 
when taking the observational results at face value as we do 
here for all datasets. After removing possible non-members 
(cf. Sec. 14.3b . a surplus of binary stars at long periods per- 
sists although less significant, which cannot be understood 
through binary-star disruption in embedded clusters. 

8. A PF of the DM type is consistent with all pre-main se- 
quence data (canonical. Lupus, UScA, Cham, UScB and 
ONC) (Hdat = 6). 

No DM type PF exists with fuM > 5 x 10 so that this 
hypothesis can be rejected with 99.95 per cent confidence. 

6. A revised IPF 

In the previous section, the tests of hypotheses |5] and |6] yield 
parameter ranges for Eqs 2] and |6] that are consistent with the 
combined data at the 5 per cent confidence level or better These 
parameter ranges are shown in Figs |5]and|6]for K-type solutions, 
and in Fig. |7] for DM-type solutions. The data sets were selected 
according to appearing least evolved, using the a-priori knowl- 
edge gained in Kl and K2. 

It is noteworthy that the K-type solutions [j] » 2.5-3,5 ~ 
40 - 70) are consistent within 5 % confidence with the Kl IPF 
(77 - 3.5,6 - 100), thus corroborating the results obtained in Kl. 
This is particularly interesting because here we used data sets in- 
cluding more than one pre-main sequence population, while in 
Kl only older Taurus- Auriga data were available but no func- 
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Fig. 7. Testing hypothesis|6]with = 1-0. Survey of a-lP,avlP 
parameter space for solutions of the DM type (Eq.|6]l to the com- 
bined canonical pre-main sequence data. Lupus and UScA pre- 
main sequence data sets. Open circles and small dots delineate 
models consistent with the data at the 1 and 5 per cent confi- 
dence level, respectively. The DM PDPF is shown as the thick 
cross. 

tional fitting to these data was performed such as here. Rather, 
the Kl IPF was obtained through inverse dynamical population 
synthesis by demanding approximate (eye-ball) fits to the pre- 
main sequence, and to the main sequence data after star-cluster 
disintegration. The present purely statistical approach confirms 
the Kl IPF (t] a 3.5,(5 ~ 100, /p^s > 0.8, with enhanced confi- 
dence probability for larger /pms)- 

A new result obtained here (Fig.|7) is that a log-normal func- 
tion in orbital period P can also be considered a parent dis- 
tribution of the same pre-main sequence data as above. These 
data are consistent, at the 5 per cent confidence level or bet- 
ter, with being drawn from a parent distribution of the DM-type, 
if this PF has parameters indistinguishable from the DM PDPF 
(alP 2.3, avlP ^ 5) but with a significantly larger binary pro- 
portion, fms = 1 -0. 

The result is thus that either the Kl IPF (Eq.|4]i, or a K-type 
solution with 6 ~ 29.52 t] (straight line in Fig.|5]l, or the Gaussian 
PF (Eq.|6]l with the above parameters can be parent distributions 
of pre-main sequence binary systems with indistinguishable con- 
fidence probabilities. This is illustrated in Fig.|8] 

7. Conclusions 

A fresh look is taken at the variations evident in the period dis- 
tribution of binary stars in various stellar populations. The most 
notable such diff'erence, noted by many authors, lies between the 
canonical pre-main sequence population of Taurus-Auriga and 
the canonical Galactic-field local G-dwarf population. The dif- 
ference is significant (SectionO, but can be understood to result 
from stimulated evolution of an initial Taurus-Auriga-type pop- 



ulation if most Galactic-field stars originate in modest embedded 
clusters (Kl). That work arrives at a possible initial period func- 
tion, the Kl IPF (Eq.EI, which is nearly flat fora >20 AU. 

The variation between pre-main sequence populations is also 
studied here, with the result that significant differences are evi- 
dent, notably between the canonical pre-main sequence popula- 
tion and the ONC and UScB samples. Stimulated evolution in 
the dense Orion Nebula Cluster efficiently depletes a Kl IPF to 
the observed data (KPM, KAH), so that the former difference 
can be readily accounted for The latter difference, which results 
from a significant apparent surplus of long-period binaries may 
indicate that the IPF of the UScB population was very different 
to that in Taurus-Auriga, as noted by Brandner & Kohler (1998). 
Noteworthy in this context is that UScA, where mostly binaries 
with small separations are present, contains several high-mass 
(B-type) stars, while no high-mass stars are present in UScB, 
which shows a preference for wider binaries. This is qualita- 
tively consistent with UScA perhaps being the result of denser 
embedded stellar groups that are already dispersed after dynam- 
ical evolution and the expulsion of residual gas by the massive 
stars. Discarding three of the widest binaries in the UScB data set 
yields a binary population consistent with the Kl IPF (Sect. 14. 31 . 
leaving the ONC as the only pre-main sequence population being 
clearly different from the canonical pre-main sequence popula- 
tion. In any case, among the data sets used in this paper the UScB 
data has the lowest number of systems surveyed for binarity, and 
the region is most likely more affected by containing dispersed 
mixed populations of stars with different distances and ages than 
the other data. Hence, the observed orbital period distribution of 
binaries in UScB remains to be confirmed. 

The (canonical) Taurus-Auriga, Lupus and UScA popula- 
tions appear least-evolved and are consistent with being drawn 
from a common parent IPF iff the binary proportion is higher 
than on the main sequence. The set of possible IPFs is shown in 
Figs.|5]|6]and|2l noting that the Kl IPF is included. Three possi- 
ble IPFs, all with indistinguishable confidence probabilities, are 
presented in Fig. [8] While no formal decision can be made based 
on !PKorDM as to which of the three IPFs plotted in Fig. |8] are to 
be preferred, the a-priori knowledge gained in Kl favours the 
Kl IPF, because it is (1) consistent with the data at the 10 per 
cent confidence level, (2) it is the precursor of the Galactic -field 
PDPF if most stars form in embedded clusters, and (3) the re- 
vised K-type IPF shown as the thin solid line in Fig. [8] leads 
to a deficit of Galactic-field binaries with IP>6 (fig. 8 in Kl). 
Inverse dynamical population synthesis will have to be applied 
to investigate if the alternative log-normal IPF (thick dashed line 
in Fig. [8]l can be made consistent with the Galactic-field PDPF 
for a reasonable library of embedded star clusters. However, by 
consulting fig. 1 of Connelley et al. (2008) it becomes readily 
apparent that a flat distribution function for a > 100 AU is pre- 
ferred over a log-normal function, as is also concluded by those 
authors and as is the case for the Kl IPF. This issue is studied in 
more detail now (Marks, Kroupa & Oh 201 1) and does suggest 
that the Kl IPF (Eq.|4| is favoured. Parker et al. (2009) come to 
a similar conclusion and stress that a > lO'* AU binaries cannot 
survive in any clustered environment. They may however form 
during cluster dissolution (Kouwenhoven et al. 2010). 

It is notable that the above three pre-main sequence popu- 
lations, and in addition the Cham and ONC populations, all are 
consistent with a monotonically decreasing PF with increasing 
IP (Figs.|2]and|8ll. This may suggest that all known pre-main se- 
quence populations have already suffered some degree of stim- 
ulated evolution, and may have begun with the same (universal) 
IPF. This notion allows reconstruction of the properties of the 
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Fig. 8. The revised IPF. In addition to the Kl IPF (thick solid curve), two possible IPF forms are displayed. These are identified as 
solutions in Figs.|5]and|7] and are shown as the thin solid line (K-type, Eq.|4]with the asymptotic solution 6 - 29.52 r];ri - 18.7 and 
/p„,, = 1.0) and the thick dashed line (DM-type, Eq.|6]with cr/P = 2.3, avZP = 5.0, /^s = 1.0). The thin dashed line is the DM PDPF 
Data points with observational error bars are as indicated in the key. These are used to test hypothesis |5]and|6]giving !P = 0.10 (thin 
solid curve), !P = 0.22 (thick dashed curve) and !P = 0.10 (thick solid curve). 



embedded clusters from which the pre-main sequence popula- 
tions could have originated (inverse dynamical population syn- 
thesis). Explicit modelling of this evolution has been performed 
for the ONC (KPM, KAH; Parker et al.2009) and sparse em- 
bedded clusters (Kl, K2, Kroupa & Bouvier 2003; Parker et al. 
2009). 

Thus, at the present it cannot be confirmed that the IPF varies 
with star-forming conditions. That is, the presently available data 
are consistent with an invariant birth or initial binary population. 
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